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We used resonant inelastic x-ray scattering to reveal the nature of magnetic interactions in 
Sr2lr04, a 5d transition-metal oxide with a spin-orbit entangled ground state and Jeff=l/2 magnetic 
momemts, referred to as 'isospins'. The magnon dispersion in Sr2lr04 is well described by an anti- 
ferromagnetic Heisenberg model with isospin one-half moments on a square lattice, which renders 
the low-energy effective physics of Sr2lr04 much akin to that in superconducting cupratcs. This 
is further supported by the observation of cxciton modes in Sr2lr04 whose dispersion is strongly 
renormalizcd by magnons, which can be understood by analogy to the hole propagation in the 
background of antiferromagnetically ordered spins in the cuprates. 



PACS numbers: 74.10.-t-v, 75.30.Ds, 78.70.Ck 

Quantum magnetism in transition-metal oxides 
(TMOs) arises from superexchange interactions of spin 

moments that depend on spin-orbital configurations in 
the ground and excited states. The array of magnetism 
in 3d TMOs are now well understood within the frame- 
work of Goodenough-Kanamori- Anderson [1], which as- 
sumes conservation of spin angular momentum in the 
virtual charge fluctuations. However, it has been re- 
cently realized that strong relativistic spin-orbit coupling 
(SOC) can drastically modify the magnetic interactions 
and yield far richer spectrum of magnetic systems be- 
yond the standard picture. Such is the case in 5d TMOs, 
in which the energy scale of SOC is on the order of 0.5 
eV (as compared to ^10 meV in 3d TMOs). For exam- 
ple, A2lr03 (A=Li,Na) is being discussed as a possible 
realization of the long-sought Kitaev model with bond- 
dependent spin interactions [2-4]. Furthermore, strong 
SOC may result in nontrivial band topology to realize 
exotic topological states of matter with broken time re- 
versal symmetry, such as a topological Mott insulator [5] , 
Wcyl semi- metal and axion insiilator [6]. Despite such 
intriguing proposals abound, the nature of magnetic in- 
teractions in systems with strong SOC remains experi- 
mentally an open question. 

In this Letter, we report on the magnetic interactions 
in a 5c? TMO Sr2lr04 with spin-orbit entangled ground 
state carrying Jou=l/2 moments [7, 8], probed by reso- 
nant inelastic x-ray scattering (RIXS) . We refer to these 
Jeff=l/2 moments as 'isospins.' On theoretical grounds, 
magnetic interactions of isospins are expected to strongly 
depend on lattice and bonding geometries. In the par- 
ticular case of corner sharing oxygen octahedra on a 
square lattice, relevant to Sr2lr04 [9] (Fig. 1(a)), it is 
predicted that the magnetic interactions of isospins are 



described by a pure Heisenberg model, barring Hund's 
coupling that contributes a weak dipolar-like anisotropy 
term [2, 10]. This is surprising considering that strong 
SOC typically results in anisotropic magnetic couplings 
that deviate from the pure Heisenberg-like spin interac- 
tion in the weak SOC limit. A compelling outcome is 
that a novel Heisenberg antiferromagnet can be realized 
in the strong SOC limit, on which a novel platform for 
high temperature superconductivity (HTSC) may be de- 
signed. 

RIXS in the last few years has become a powerful tool 

to study magnetic excitations [11]. We report measure- 
ment of single magnons using hard x-rays, which has 
complementary advantages over soft x-rays as detailed 
later on. The RIXS measurements were performed us- 
ing the spectrometer at the 9-lD and the MERIX spec- 
trometer at the 30-ID beamline of the Advanced Photon 
Source. A horizontal scattering geometry was used with 
the TT incident photon polarization. A spherical diced 
Si (844) analyzer was used. The overall energy and mo- 
mentum resolution of the RIXS spectrometer at the Ir 
L3 edge (R:!ll.2 keV) was ^130 meV and ±0.032 l"!, 
respecitvely. 

Sr2lr04 has a canted antiferromagnetic (AF) struc- 
ture [8] with Tn « 240 K [12], as shown in Fig. 1(b). 
Although the 'internal' structure of a single magnetic 
moment in Sr2lr04, composed of orbital and spin, is 
drastically different from that of pure spins in La2Cu04, 
a parent insulator for cuprate superconductors, the two 
compounds share apparently similar magnetic structure. 
Figures 2(a) and (b) show the dispersion and intensity, 
respectively, of the single magnon extracted by fitting 
energy distribution curves shown in Fig. 3(a). We high- 
light three important observations. First, not only the 
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FIG. 1. (a) Due to a staggered in-layer rotation of oxygen oc- 
tahedra, Sr2lr04 has four Ir02 layers in the unit cell [9], which 
coincides with the magnetic unit cell, (b) Isospin one-half 
(Jeff — 1/2) moments lie and are canted in the Ir02 plane [8]. 



dispersion but also the momentum dependence of the 
intensity show striking similarities to those observed in 
the cuprates by inelastic neutron scattering, for instance 
in La2Cu04 [13]. This provides confidence that the ob- 
served mode is indeed a single magnon excitation [14- 
17]. Using hard x-ray RIXS allows mapping of an entire 
Brillouin zone within only a few degrees of 90° scatter- 
ing geometry so that the spectrum reveals the intrinsic 
dynamical structural factor with minimal RIXS matrix 
element effects. Second, the measured magnon disper- 
sion relation strongly supports the theories predicting 
that the superexhange interactions of isospins on a square 
lattice with corner-sharing octahedra are governed by a 
SU(2) invariant Hamiltonian with AF coupling [2, 10]. 
Third, the magnon mode in Sr2lr04 has a bandwidth 
of ~200 meV as compared to ^300 meV in La2Cu04 
[13] and Sr2Cu02Cl2 [18], which is consistent with en- 
ergy scales of hopping t and on-site Coulomb energy U 
in Sr2lr04 being smaller by roughly 50% than those re- 
ported for the cuprates [10, 19, 20]. 

For a quantitative description, we have fitted the 
magnon dispersion using a phenomenological J-J'-J" 
model [22]. Here, the J, J', and J" correspond to the 
first, second, and third nearest neighbors, respectively. 
In this model, the downward dispersion along the mag- 
netic Brillouin zone from (7r,0) to (7r/2, 7r/2) is accounted 
for by a ferromagnetic J' [13, 22]. We find J=60, J'=- 
20, and J"=15 meV. The nearest neighbor interaction J 
is again smaller than found in cuprates by roughly 50%. 
The fit can be improved by including higher order terms 
from longer-range interactions, which were also found to 
be important in the case of Sr2Cu02Cl2 [18]. However, 
here we do not pursue this path because, as we show be- 
low, another kind of magnetic modes in Sr2lr04, which 
is not present in cuprates, may affect the magnon disper- 
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i^lG. 2. (a) Blue dots with error bars show the smgle magnon 
dispersion extracted by fitting the energy loss curves shown 
in Fig. 3(b) [21]. The magnons disperse up to ~205 meV at 
(7r,0) and 110 meV at (7r/2,7r/2). Solid purple line shows the 
best fit to the data with J=60, J' =-20, J" = 15 meV. (b) Mo- 
mentum dependence of the intensities showing diverging in- 
tensity at (tt, tt) and vanishing intensity at (0,0). Inset shows 
the Brillouin zone of the undistorted tetragonal (I4/mmm) 
unit cell (black square) and the magnetic cell (blue square), 
and the notation follows the convention for the tetragonal unit 
cell, as, for instance, in La2Cu04 . 

sion. 

Ultimately, characterizing the magnon mode is impor- 
tant because it strongly renormalizes the dispersion of a 
doped hole/electron, and is believed to provide a pair- 
ing mechanism for HSTC. We now show that Sr2lr04 
supports an exciton mode, which gives access to the dy- 
namics of a particle propagating in the background of 
AF ordered moments even in an undoped case. Fig- 
ures 3(a) and (b) show the energy loss spectra along 
high symmetry directions. No corrections to the raw 
data such as normalization or subtraction of the elas- 
tic contaminations have been made. Another virtue of 
using hard x-ray is that by working in the vicinity of 
90° scattering geomeotry elastic (Thompson) scattering 
can be strongly suppressed. In addition to low-energy 
magnon branch (< 0.2 eV), we observe high-energy ex- 
citations with strong momentum dependence in the en- 
ergy range of 0.4~0.8 eV. This mode is superimposed on 
top of a continuum generated by particle-hole excitations 
across the Mott gap (estimated to be «0.4 eV from op- 
tical spectroscopy [23]). This is schematically shown in 
Fig. 3(c). Taking the second derivative of the raw data 
de-emphasizes the intensity arising from the particle-hole 
continuum and reveals a clear dispersive feature above 
0.4 eV, as shown in Fig. 4(a). The energy scale of this 
excitation coincides with the known energy scale of spin- 
orbit coupling in Sr2lr04 (Cso ^^0.5 eV) [7], and thus 
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FIG. 3. (a) Energy loss spectra recorded at T=15 K, well below the Tn «240 K [8, 12], along a path in the constant L=34 
plane. The path was chosen to avoid the magnetic Bragg peaks, which appear at two of the four corners of the unfolded unit 
cell (black square) shown in the inset (where the same conventions as in Fig. 2 are used), (b) Image plot of the data shown in 
(a), (c) Schematic of the three representative features in the data. 



we assign it to intra-site excitations of a hole across the 
spin-orbit split levels in the t2g manifold, i.e. from the 
Jgff=l/2 level to the one of the Joff=3/2 levels [7, 14, 21] 
(see Fig. 4(c)). We refer to such an excitation as 'spin- 
orbit exciton' [24]. 

The dispersion of spin-orbit exciton with a bandwidth 
of at least 0.3 eV implies that this local excitation can 
propagate coherently through the lattice. Our model of 
the spin-orbit exciton starts from a recognition that the 
hopping process is formally analogous to the problem of 
a hole propagating in the background of AF ordered mo- 
ments, which has been extensively studied in the context 
of cuprate HTSC [25]. Although the spin-orbit exciton 
does not carry a charge, its hopping creates a trail of 
misaligned spins and thus is subject to the same kind 
of renormalization by magnons as that experienced by a 
doped hole [26] . It is well known that the dispersion of a 
doped hole in cuprates has a minimum at (7r/2,7r/2) [27], 
i.e., at the AF magnetic Brillouin zone boundary. Since 
Sr2lr04 has a similar magnetic order [8], it can be un- 
derstood by analogy that the dispersion of the spin-orbit 
exciton should also have its minimum at (7r/2,7r/2). 

The overall bandwidth is determined by the parame- 
ters involved in the hopping process, which is depicted 
in Fig. 4(c) in the hole picture. It involves moving an 
excited hole to a neighboring site, which happens in two 
steps. First, the excited hole in site i hops to a neighbor- 
ing site j, generating an intermediate state with energy 
U' , which is the Coulomb repulsion between two holes at 
a site in two different spin-orbital quantum levels. Then, 
the other hole in site j hops back to site i. Thus, the en- 
ergy scale of the dispersion is set by 2tAAtBB/U' , which is 
of the order of the magnetic exchange couplings. In fact, 
these processes lead to the superexchange interactions re- 
sponsible for the magnetic ordering, but here they involve 
both the ground state and excited states of Ir ions. 



Technically, the spin-orbit exciton hopping can be de- 
scribed by the following Hamiltonian 

where i indicates the lattice site, b (b^) is the magnon 
annihilation (creation) operator, and X denotes the spin- 
orbit exciton that carries a quantum number a belonging 
either to the B or C doublet in the Joff=3/2 manifold 
(see Fig. 3(c)). From this expression, the analogy with 
the case of a moving hole is apparent; in place of the 
hopping t for the doped hole, we have an effective spin- 
orbit hopping matrix W with its overall energy scale set 
by = 2t^/U. 

We calculated the spin-orbit exciton dispersion by eval- 
uating its self-energy matrix expressed as 

where z is the coordination number and M denotes the 
vertex [21], using the actual experimental magnon dis- 
persion relation for Wq as shown in Fig. 2(a). The only 
adjustable parameter is W, which only contributes to 
the overall scaling of the dispersion. The model correctly 
captures the main features of the data: the locations of 
extrema in the dispersion (Fig. 4(a)), nearly momentum 
independent integrated spectral weight (Fig. 4(b)), and 
the intensity relative to the magnon intensity (Fig. 4(b)). 
The theory described above is developed in Supplemen- 
tary Material. 

Our measurement of the spin-orbit exciton dispersion 
has important implications in modeling 5ci transition- 
metal oxides with strong SOC. First, it shows that 
not only the Jeff =1/2 states are localized but also the 
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FIG. 4. (a) Second derivative of the data shown in Fig. 3(b), 
overlaid with the calculation of the dispersion (red solid and 
dashed lines). The upper branch (dashed line) has less spec- 
tral weight than the lower branch (solid line) [21]. W=63 
meV was chosen in the calculation, which is within the es- 
timated range [21]. (b) Comparison of the experimental 
and theoretical integrated spectral weight for the two spin- 
orbit exciton modes normalized to the single magnon mode 
at (7r,0). (c) Schematic of the spin-orbit exciton hopping in 
the hole representation. By the RIXS process, the hole in the 
i site is excited to the B/C doublets, which are the Jefi=3/2 
quartet split by the tetragonal crystal field. This excited hole 
hops to the neighboring site j with the intermediate energy 
of U . The other hole in the j site hops back to the i site, 
thereby completing the spin-orbit exciton hopping processes 
and also creating a magnon (blue wavy line). 



Joff=3/2 states largely retain their atomic-like charac- 
ter. In a contrasting model, in which Joff=3/2 states 
form an itinerant band and only the Joff=l/2 states are 
localized, much akin to the orbital-selective Mott tran- 
sition scenario [28], one expects to see only an electron- 
hole continuum resulting from the independent propa- 
gations of a hole and an electron. Instead, we see the 
spin-orbit exciton coexisting with the particle-hole con- 
tinuum; a duality of atomic and band nature of the same 
hd electrons. Second, the existence of the virtually bound 
J(,ff=3/2 states only ~0.5 eV above the ground state im- 
plies that the superexchange interactions entail multior- 
bital contributions. Thus, even for an apparently single 
orbital Jch— 1/2 system such as Sr2lr04, the magnetic in- 
teractions are multiorbital in character, a fact that must 
to be taken into account in any quantitative model. 

Despite such important differences in the high energy 
scale, our measurement of the magnon spectrum high- 
lights the similarities with cuprates in the low energy 
effective physics - a rare realization of (iso)spin one-half 
moments on a square lattice with Heisenberg SU(2) in- 
variant interactions and comparable magnon bandwidth. 



Further, from the observed spin-orbit exciton dispersion, 
we may expect that a doped hole/electron in Sr2lr04 
will display the same dynamics as that observed for a 
doped hole/electron in the cuprates. The phase diagram 
of lightly doped Sr2lr04 has just begun to be revealed 
experimentally [29, 30]. Although superconductivity has 
not yet been reported, some anomalies that bear strong 
resemblance to cuprates such as T-linear resistivity have 
been seen [30, 31]. Only further study will tell if doping 
can drive Sr2lr04 superconducting. 
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I. RIXS PROCESS AND EFFECTIVE 
OPERATORS 

We describe the RIXS process and derive the effective 
RIXS operators for the particular case of Sr2lr04^^. The 
states we use are A± for Joff = 1/2 ground state doublet 
(where the sign takes the role of a isospin degree of free- 
dom) and B± and C± Kramers doublets for Jeff = 3/2 
excited states (see Fig. SI). The states A±, B±, and C± 
depend on spin orbit coupling A and tetragonal crystal 
field A^^. The energy splitting A between the A± states 
and the B±/C± subsystem has a dominant contribution 
from spin-orbit coupling |A but also may be contributed 
by coupling to the lattice; The charge density distribution 
in the excited Joff = 3/2 states is different from that in 
the ground state, and hence the orbital flip may cost ad- 
ditional energy. In terms of electron counting, the B and 
C doublets are completely occupied, while the A level is 
half filled. We now switch to the hole notation, in which 
B and C are empty. In the A subsystem, superexchange 
interactions stabilize an alternating magnetic (isospin) 
order. We consider the low-energy Hilbert space with 
singly occupied sites with an empty B/C subsystem, re- 
ferred to the hole notation. In the Mott insulating ground 
state, all these levels are coupled by virtual hoppings re- 
sulting in a multiorbital superexchange interaction oper- 
ating within the Hilbert space spanned by Jeff =1/2 and 
JefT = 3/2 states. 

Figure SI illustrates the RIXS process in the hole rep- 
resentation. The hole in the A level is excited to the 2p 
core levels, also consisting of three doubles (a, 5, and c), 
by absorbing a photon. Subsequently, the excited inter- 
mediate state decays back to the t2g manifold, emitting 
a photon. The excited hole decaying back to the A level 
with opposite a corresponds to an effective isospin flip, 
and the B/C level to an effective spin/orbital flip rela- 



Cround state Intermediate state 



Final Excited state 



L2 (2p,/2) a 
L3 {2p„2) b/c 



b/c =(= 



photon in 



photon out 



J,ff.3/2 B/C = : : . 

A _(_ + 



J,ff=l/2 



B/C 
A 



b/c === 

B/C =^= 

A 



FIG. SI. The RIXS process in the hole representation. 



five to the orbital/spin. To describe the RIXS process 
mathematically, we define operators Xa-, and Z^- for 
annihilating a hole in the yz, zx, and xy orbital, respec- 
tively. These operators can be expressed in terms of the 
A^, Ben ^nd Co- Kramers doublets as 



V2 



(cos 6A_„ - sin 9B_„ + C_ 



Y„ = — (cosM_o - sin6'S_o - C- 
Z„ = sin 6A„ + cos 9B„ ■ 



(1) 

(2) 
(3) 



Likewise, we define operators for creating a hole in the 
2p manifold as 
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When the spin-orbit coupling is strong and dominates 
the crystal field effects, as realized in Sr2lr04^^), then 
cos (j) = cos 9 — 

The effective RIXS amplitude A^ can be defined as 



A^<^{!\Y^e^^■^D\l)D{e)\^), 



(7) 



R 



where D given by 

B = ^Y^{tyZ^ +e^y'^)„X„^((i^x'^ +t,^z^)„Y„+(t.j,y'^ +<iyx'^)„Z„ 

(8) 

is the dipole transition operator for exciting the valence 
hole into the p core level, a function of incoming photon 
polarization e. Likewise describes the decay and is a 
function of the outgoing photon polarization e'. 

Integrating out the p manifold and using Eqs.1-3, the 
effective RIXS operator R = {£')D{e) is obtained. 
This operator annihilates a particle in the A doublet and 
creates a particle in the A doublet or in the B/C subsys- 
tem: 

^ = E + l3„,a'Bl + l.,.,Cl) A, 



Ra + Rb + Rc 



(9) 



The coefficients ar,a' , Paa' , and jaa' depend on the scat- 
tering geometry and determine the intensities. For the 
L2 edge, they all vanish in the limit of cubic symme- 
try; shown below are the results for the L3 edge relevant 
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here. The first term i?^ describes processes involving 
the A doublet only and is thus sensitive to magnons. Ex- 
pressing A^A in terms of isospin operators, Ra can be 
recast into the form 



Ra 



(10) 



where P=e' x e and S denotes the isospin within the A 
doublet. The minus sign in front of the component 
reflects the fact that the effective g factor of this doublet 
is opposite to that of real spin because of large orbital 
moment^^. In principle, the Ra can be transformed to 
a scalar product form, by a simultaneous sign change of 
the S^,Sy components. 

Performing similar calculations, the matrix elements 
for the processes described by -Rb and Rc, which create a 
spin-orbit exciton by promoting a hole from the A states 
to the B/C levels, can be derived as 



1 



<^{Ty-lPy)+i{Tx + lPx) 



(11) 
(12) 

(13) 



2V2 

-^{Q2 + i(TTz), 
2^ [-a(T, + Py) + i{T, - P,)] , (14) 



with 



Q2 
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— £y£y - 
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Q3 








• e — > 


Tx 


= e'yCz + 
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<^'z^x, 


Tz 


= s^^ + 


1 

^x^V 



(15) 
(16) 

(17) 

(18) 
(19) 



II. MODELING THE SPIN-ORBITON 
DISPERSION 

In a leading approximation (neglecting the intensity 
transfer between the exciton and magnons due to a mode- 
mode coupling), the combined spectral weight oiRs+Rc 
transitions averaged over the outgoing polarization e' is q- 
independent, and is given by Ea,<r'(l/^<^,<^'l^ + h<y,a'\'^) = 
|. The mode- mode coupling is expected to bring only 
weak q dependence to the exciton intensity since the re- 
lated corrections (of the order of J/ A) are small. This 
agrees well with the experimental fact that the spectral 
weight of the high-energy mode is nearly momentum in- 
dependent. 

On the other hand, the magnon intensity which is de- 
termined by a dynamical spin structure factor strongly 

depends on q, see Eq. 10. For the present experimen- 
tal setup, the RIXS signal is determined by a following 
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FIG. S2. Spectral weight of the magnon excitation divided 
by the weight of the excitonic sector at higher energy from 
the calculation. 



combination of the transverse xx, yy and longitudinal zz 
components of the isospin susceptibility: 



97r^2^ 



(20) 



from which the intensities of single-magnon I^^ (q) and 

two-magnon I^^ (q) processes can be calculated. (In this 
expression, we neglected weak q dependence of the po- 
larization factors). Within a linear spin- wave theory, we 
find: 




Aa+B' 



^ / ^k^k-l-q — ^k^k-l-q _ ^ 
1, V 



Wk^k-l-q 



(21) 



with Aq, Bq and magnon dispersion Wq given below in 
Eqs. 29 and 30. 

The relative intensity of the magnon and excitonic 
modes is a parameter-free result and provides an impor- 
tant test for the consistency of the model and the assign- 
ment of the observed modes. The ratio of the intensities 
of the magnon and the excitonic high-energy modes is 
shown in Fig. S2. For q = (7r,0), for instance, theory 
predicts the ratio ~ 0.3 which compares well with the ex- 
perimental value 0.26 of the relative intensities of the low- 
and high-energy modes. We now discuss the origin of the 
dispersive features in the excitonic channel. The exciton 
created by Rb and Rc delocalizcs due to virtual hop- 
pings of electrons as depicted in Fig. 3c. In essence, this 
propagation is due to a multiorbital superexchange inter- 
action which is responsible for magnetic ordering of the 
A Kramers doublets, but involving here both the ground 
>/cff=l/2 state and excited Jeff— 3/2 states of Ir ions. As 
can be seen in Fig. 3c, the two holes (drawn as isospin ar- 
rows) are always in different orbitals. (Virtual processes 
involving doubly occupied orbitals exist, but do not con- 
tribute to exciton delocalization and hence will not be 
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considered here.) If we neglect the Hund's rule couplmg, 
this implies that the spin of the excited electron plays no 
role, i.e., the hole-excitation pair can be seen as just a 
spinless hole, as in the "usual" t-J model derived from 
the one-band Hubbard model. The hole has, however, an 
orbital flavor: B or C. 

In the two-dimensional space spanned by these two or- 
bitals, the effective hopping matrix for the exciton reads 



x/y 



'2tAA 
U' 



±tBC 

tec 



(22) 



where the inter-orbital element changes sign between the 
X and y directions. One clearly sees the similarity to the 

.2 

superexchange: the effective exciton hopping is oa 
i.e., of the order of J instead of t. [The factor 2 arises 
because the second step in Fig. 3c can also occur first.] 

While the spin of the excited B/C state does not en- 
ter the Hamiltonian, the exciton creates a hole in the 
ordered isospin background of the A system below the 
magnetic ordering temperature. Following Ref. S5, this 
hopping connected to magnon creation/annihilation can 
be written as 



H]<\n — 



E 



(23) 



where we denoted a spin-orbit exciton by Xi a- It carries 
an index i indicating a lattice site and a quantum number 
a indicating whether it is in the B or C doublets. We 
see that a hopping of X particle is accompanied by a 
creation or annihilation b of magnon excitations. W"j^ 
is the orbital- and direction-dependent effective hopping 
amplitude, see Eq. 22. We consider here only the nearest- 
neighbor hoppings. 

Comparison of our effective Hamiltonian to that of a 
hole moving in a quantum antiferromagnet reveals that 
the two problems are equivalent, even though we do not 
describe a real charge moving at an energy scale of t, 
but exciton motion with a strongly renormalized hopping 
matrix W"'^ = Wt'^'^. The energy scale W = 2t\j^/U' 
quantifies the exciton hopping amplitude, and Eq. 22 be- 
comes 



w: 



a/3 

■/y 



Hsc/tAA 



HBc/tAA 
tcc/tAA 



Q/3 



(24) 



where t°'^ are numbers of the order of one {tsB/tAA — 

j,tBc/tAA = tcc/tAA = I in the cubic limit 

used here; the ± sign depends on x/y bond direc- 
tion). For nearest-neighbor hopping of the Jeff ==1/2 
states, tAA — 0.26 eV is suggested from band structure 
calculations^^" and using a representative value of in- 
terorbital Coulomb interaction U' ^ 2 eV, one estimates 
1^- 60- 70 meV. 

Each hopping process creates or destroys magnons. We 
follow Ref. S5 and treat the AF background in the lin- 
ear spin-wave theory. After the Bogoliubov transforma- 
tion, the correlated hopping of exciton (accompanied by 



(ji/2,n/2) (jt,0) 



FIG. S3. Dispersion of the hole-excitation pair. Color shading 
shows the spectrum convoluted with a Lorentzian broadening 
with width 5=100 meV and weighted with their intensity, see 
Eq. 32, the lines show the eigenvalues of A-|-Ek, see Eq. 31. 
Parameters: W=63 meV, Ji=60 meV, J2=-20 meV, J2=15 
meV, and A =0.78 eV. 



a magnon creation/annihilation) reads as follows: 



k,q,Q,/3 



(25) 



Except for the orbital index, this is analogous to the t-J 
model for a doped hole motion in AF-ordered cuprates. 
Here, X^^. {X^^J creates (annihilates) an exciton with 
momentum k and orbital flavor a, while (oj^) creates 
(annihilates) a magnon, and z = 4 is the coordination 
number. A momentum- and orbital-dependent vertex has 
been introduced as 



a/3 
k,q 



I a/3|/ ap , Oit 



(26) 



absorbing the sign change of t bc into the following form- 
factors: 



7k 



a0 



1 f {cos + COS ky) {cos — cos ky) 



2 I (cos kx — cos ky) (cos k^ + cos ky) 



(27) 



Wq and Wq are Bogoliubov factors and given by the rela- 
tions 




sgn^q /Ac 
\ 0J^ 



- 1 



(28) 



with 



Aq = 2( Ji - J2 - J3 + J2 cos Qx COS qy) + J3(cos + cos 2qy), 



Br, 



Ji{cosqx + cosqy), 



and the magnon energy 



.q = - Bl 



(29) 



(30) 



The parameters Ji, J2 and J3 denote nearest-neighbor, 
next-nearest neighbor, and third-neighbor couplings of 
the isospins Sa- 
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The total Hamiltonian describing the exciton motion 
in the AF-ordered background contains the kinetic en- 
ergy ifkin as well as the energy of the magnon sector 
Hi^ag — J2q ^q^q^q- Sincc W'^y^ is of the order of the J's 
rather than the bare hopping t, the present case is slightly 
different from the "usual" t-J model obtained as a large- 
U limit of the Hubbard model, where t ^ J. The latter 
describes a strong coupling of holes to magnons, while the 
present case has a coupling comparable to the magnon 
energy, ~ J, since both parameters originate from the 
same multiorbital superexchange interaction. We conse- 
quently use perturbation theory to describe the exciton 
moving in the AF background. Second order perturba- 
tion theory leads to an energy correction 

i:lP = -{zWfY. "^'^ '^'^ (31) 

with the vertex M^^ defined in Eq. 26. This sum can 
be evaluated numerically, and with parameters W — 
63 meV and A = 0.78 eV we find the dispersive feature 
shown in Fig. S3. 

In Fig. S3, the lines give the two eigenvalues of S^*^, 
which determine the excitation energies. The weight of 
the poles can be obtained from the prefactors given in 
Eq. 11 and the eigenvectors cos 0i(2)k7 sin0]^(-2)k corre- 
sponding to the two eigenvalues of The spectral 
weight of each pole (labeled by 1 and 2) is given by 

/l(2) (k) = [\Pa,a'f COS^ 6*1(2) + \la^a' f siu^ 6*1(2) 

a. a' 

(32) 

We recall that while the partial intensities /i(k) and 
I2 (k) are momentum dependent (via the cos and sin 0k 
functions) affecting thereby the energy/momentum in- 
tensity map, the total spectral weight Ji -I- I2 is nearly 
constant (= |) over an entire Brillouin zone. 

Corrections to the dispersion curves can be calculated 
within a standard self-consistent Born approximation 
(SCBA). However, the basic features (such as location of 
minima at the magnetic Brillouin zone boundary, max- 
ima at r and magnetic Bragg points, overall bandwidth 
of the order of magnon bandwidth, exciton/magnon rel- 
ative intensities) are expected to remain unchanged since 
these observations are intrinsic to the model. Based on 
a well-studied case of hole motion in quantum antiferro- 
magnets, we may expect that the major effect of SCBA 
calculations will be in a substantial broadening of the 
exciton lines. 



III. EXTRACTION OF THE SINGLE MAGNON 
DISPERSION 

Figure S4 shows the low-energy part of the spectra 
fit with three peaks along the magnetic Brillouin zone. 




(ji,0) 



(ji/2,3t/2) 
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FIG. S4. The single magnon mode extracted from peak fit- 
ting. The square markers show the experimental data and 
the red solid lines shows the fit to the data. The red shaded 
peak corresponds to the single magnon peak. The dashed 
black and blue fine correspond to the elastic and two magnon 
peaks, respectively. 

The three peaks account for the elastic, single (shaded 
peaks), and double magnon peaks. The elastic and sin- 
gle magnon peaks were constrained to have a Gaussian 
width of 130 meV, which is the instrumental resolution 
determined from the elastic peak at a nonresonant condi- 
tion. The high-energy modes above «0.4 eV were also fit 
with Gaussian functions to account for the tail extend- 
ing to the magnon modes. The double magnon peaks are 
not clearly resolved, especially near (7r,0), which results 
in some uncertainties in the single magnon peak posi- 
tions and amplitudes. However, regardless of the fitting 
procedure used the single magnon peak position was not 
affected by more than 5 meV. 
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